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Abstract. 

Positronium is an ideal system for the research of the bound state QED. New precise measure- 
ment of orthopositronium decay rate has been performed with an accuracy of 150 ppm, and the 
result combined with the last three is T.CWOliO.OOOVjtu -1 . It is the first result to validate the 2nd 
order correction. The Hyper Fine Splitting of positronium is sensitive to the higher order corrections 
of the QED prediction and also to the new physics beyond Standard Model via the quantum oscilla- 
tion into virtual photon. The discrepancy of 3.5(7 is found recently between the measured values and 
the QED prediction (0(a 3 )). It might be due to the contribution of the new physics or the systematic 
problems in the previous measurements: (non-thermalized Ps and non-uniformity of the magnetic 
field). We propose new methods to measure HFS precisely without the these uncertainties. 
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INTRODUCTION 

Positronium (Ps), the bound state of an electron and a positron, is a purely leptonic 
system and the triplet (l 3 Si) state of Ps, orthopositronium(o-Ps), decays slowly into 
three photons. It is good sample to measure decay rate directly and precisely. The o- 
Ps has an energy higher than the single state (l 1 ^) °f P s > parapositronium(p-Ps). The 
difference of the energy level is called as Hyper Fine Splitting (HFS) and is significantly 
larger( about 203 GHz) than the hydrogen-atom (1.4 GHz). Since o-Ps has the same 
quantum number of photon, the quantum oscillation o-Ps — > y* — > o-Ps contributes the 
HFS (87 GHz), and the oscillation has the good sensitivity to the new physics beyond 
Standard Model. 

Precise measurements of the decay rate of o-Ps and the HFS give us direct information 
about quantum electrodynamics(QED) in bound state. The higher order calculations in 
the bound state have be developed recently, and the we have a good chance to compare 
the experimental results with the QED predictions directly. 

There are two topics in this note. We report the first test of the 2nd order calculation 
of the orthopositronium decay rate in the first part. We perform the accurate measure- 
ment (error <150ppm) of the o-Ps decay rate, and the 2nd order calculation is directly 
validated for the first time. In the second part, we propose the new methods to measure 



the HFS. There is discrepancy(3.5<7) between the prediction(0(oc 3 )) and the precise 
measured values. We point out the possible two systematic errors in the previous mea- 
surements and propose the new methods without these systematic uncertainties. 

ORTHOPOSITRONIUM DECAY RATE 

History and current status 

Three precise measurements 01 0,11] of the o-PS decay rate were performed, in which 
reported decay rate values much larger, i.e., 5.2 - 9.1 experimental standard deviations, 
than a QED prediction 0] (7.039934(10) Hs~ l ) corrected up to 0(a 2 ). This discrepancy 
has been referred as 'orthopositronium lifetime puzzle', and was long-standing problem. 

As some fraction of o-Ps inevitably results in 'pick-off annihilation due to collisions 
with atomic electrons of the target material, the observed o-Ps decay rate is a sum 
of the intrinsic o-Ps decay rate A -p s and the pick-off annihilation rate into 2y's, A p;c £, 
i.e., 

hobs = ^3y + hpick- ( 1 ) 

In the old measurements 01 0, 0], A ^'s were measured by varying the densities of 
the target materials, size of the cavities and also the entrance aperture of the cavities. 
The extrapolation to zero density or aperture was expected to yield the decay rate in a 
vacuum, X^y, under the assumption of the constant o-Ps velocity. We pointed out[j5|, @] 
that this assumption was the common/serious systematic uncertainties in these results, 
since it takes much time that Ps is well thermalized. 

We have proposed the following entirely new method^, which is free from above- 
mentioned systematic error. The energy distribution of photons emitted from the 3-body 
decay is continuous below the steep edge at 51 1 keV, whereas the pick-off annihilation is 
2-body which produces a 5 1 1 keV monochromatic peak. Energy and timing information 
are simultaneously measured with high-energy resolution germanium detectors such that 
Xpi C k{t) I Xiy can be determined from the energy spectrum of the emitted photon. Once 
a precise thermalization function is obtained, X p i c k{t) will contain all information about 
the process. The population of o-Ps at time t, N(t) can be expressed as 

N(t) =^exp (-A 3y jf (l + i*^) dt'^j . (2) 

Providing the ratio is determined as a function of time, the intrinsic decay rate of o-Ps, 
Xt,j, can be directly obtained by fitting the observed time spectrum. 

We obtained decay rate of 7.0398(29) /i[6], are consistent with the first order QED 
calculation, and quite differ from the previous results as shown in Fig. 1. This result 
are confirmed by two different and more accurate (200ppm) measurements [0. Si and the 
'lifetime puzzle' has been solved. Now the interesting is focused on the validation of the 
higher order correction itself. 

Non relativistic QED approximation has been developed recently, and it is useful 
to calculate higher oder correction of the bound state. The second order correction [0], 
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FIGURE 1. Historical plots of the measured o-Ps decay rate 



whose contribution is about 160ppm, has been performed in 2002, and the more accurate 
measurement is necessary to examine the 2nd order correction. 



Experimental setup 

Figure 2 shows a diagram of the experimental setup!!; A 68 Ge — Ga positron source 
(dia., 10 mm) with the strength of 0.3/iCz, being sandwiched between two sheets of 
plastic scintillators(NE102 thickness=200/im) and held by a cone made of aluminized 
mylar. The cone was situated at the center of a cylindrical vacuum container made 
of 1 mm-thick plastic scintillators and glass, being filled with Si02 aerogel (RUN-I) 
or powder (RUN-II), and evacuated down to 1 x 10~ 2 Torr. Density of the aerogel 
and powder is 0.03g/cc both, and the surface of the primary grain are replaced into 
hydrophobic in oder to remove the electric dipole of the OH-. The sizes of the primary 
grain are 10 and 16 nm for the aerogel and powder, respectively. Thus the mean free path 
of the collision between positronium and S1O2 are different, and it makes difference in 
the thermalization process and the pickoff probability. 

The emitted positrons pass through thin scintillator, transmitting a light pulse to a 
photomultiplier(Trigger PMT) and forming Ps when stopped in the silica powder or 
aerogel. Since the kinetic energy of the positron is high (E ene j=l.9MeV: 68 Ge — Ga) 
about 40% of the positron pass through the silica and annihilate at the vacuum container. 
If the positron deposits energy(>100 keV) in the 1 mm-thick plastic scintillators, inside 
wall of the vacuum container, the produced trigger signal is vetoed to remove these 
annihilate events in the wall. 

Three high-purity coaxial germanium detectors (Ortec GEM38195) precisely mea- 
sured the thermalization process. Energy resolutions were measured using several line 
7 sources, with typical resultant values of 0.5 keV at 514 keV. Four YAP(YA1C>3: Ce 
doped: size is 50x50x33 mm) scintillators simultaneously measured the time and en- 



Alnlnind mflar <t-»*a) 




Tr iBaor 1 2 inch PUT 



PlHllo MrlnlH Inter 
for E r buaf 

(i-ZQOi.*, H-ICt**) Anil Ir,fia«r 1 inchMWf 



I iKi.r culd* 

^IbIId acini 1 1 tator 
Far rtpii I - k 1 1 u«r 



TAP M I rtl M lotor 
t Mum x 5flnm * Hnm ) 



FIGURE 2. Schematic diagram of apparatus.! 



ergy information from each decay. The good time resolution of (1 nsec) is obtained for 
Ey > 150 keV. Since the decay constant of the YAP is fast(30 nsec) and the slow compo- 
nents of the scintillation is negligible small[9], the event pileup probability is small and 
the time walk correction can be performed perfectly as function of the deposited energy. 



Analysis 

The ratio ^pickify / ^iy lS determined using the energy spectrum measured by the 
Ge detectors. The energy spectrum of the o-Ps decay sample, referred to as the o- 
Ps spectrum, is obtained by subtracting accidental contributions from the measured 
spectrum. The 3y-decay continuum spectrum is estimated using Monte Carlo simulation 
in which the geometry and various material distributions are reproduced in detail|9|]. For 
every simulated event, three photons are generated according to an order- a-corrected 
energy spectrum. Successive photoelectric, Compton, or Rayleigh scattering interactions 
of every photon are then followed through the materials until all photon energy is 
either deposited or escapes from the detectors. The response function of the detectors 
is determined based on the measured spectrum of monochromatic y-rays emitted from 
152 Eu, 85 Sr, and 137 Cs, with this function being used in the simulation [@]. These material 
and detector effects are taken into account precisely, and the 3y- spectrum is obtained. 
These are continuous distributions, and is normalized to the observed o-Ps spectrum 
with the ratio of event numbers within the region (480-505 keV). 

Figure 3(a) shows the o-Ps spectrum and 3y spectrum observed in RUN-II (powder). 
The pick-off 2y spectrum, which is obtained from the o-Ps spectrum after subtracting the 
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FIGURE 3. (a) Energy spectrum of o-Ps decay y's obtained by Ge detectors. The solid line represent 
data points in a time window of 60 — 700 nsec, and the dotted line shows the 3 /-decay spectrum calculated 
by the Monte Carlo simulation. Pick-off spectrum obtained after subtracting the 3 7 contribution from the 
o-Ps spectrum is superimposed, (b) The ratio X p i C k{t) / Xj,y are plotted as a function of time for RUN- 
Il(powder). Only the statistical errors are shown and the solid lines represent best fit results obtained. 



37-spectrum, is also superimposed in the same figure. The Xpick/^y can be calculated 
directly from the ratio of the pick-off 2y and 3 7 spectrum, the ratio does not depend 
on the absolute values of the detection efficiencies. The calculations of the ^pick/^y 
are performed for the various time windows, and its time dependence is shown in 
Fig. 3(b). Horizontal axis of the figure is the time between the positrons emission and 
decay, this slop show the thermalization process of the positronium. It takes much 
time than the lifetime to well thermalized, as we have already shown in the previous 
measurements [0, 0] • 

Figure 4 shows the time spectrum of the YAP scintillators with an energy window 
above 150 keV, in which a sharp peak of the prompt annihilation is followed by the 
exponential decay curve of o-Ps and then the constant accidental spectrum. Time reso- 
lution is 1.2 nsec is obtained and the o-Ps curve is widely observed over 1.0 fls. We fit 
resultant time spectrum using the least square method, i.e., 

(3) 

where £ p i c k and £3^ are the detection efficiencies for pick-off annihilation and 3y decays, 
and R st0 p is an experimental random counting rate representing the fact that time interval 
measurement always accept the first 7 as a stop signal. ^ P i c k/^3y is about 1% due to the 
low-density of the S1O2 powder or aerogel, i.e., the ratio of error propagation to decay 
rate is suppressed by a factor of 100. 

The obtained decay ratesH are A 3y = 7.03876 ±0.0009 (stat.) [is' 1 for RUN-I and 
7.04136 ±0.0009 (stat.) /ls~ l for RUN-II, which are consistent with each other. 

Estimates of various systematic errors are summarized in Table 1 and the details 
are shown in reference[9]. We just mention here about errors related to simulation and 
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FIGURE 4. Timing spectrum measured by the YAP scintillators within an energy window above 
150 keV for RUN-II. The figure in the right corner shows the magnified view of the prompt peak. 



TABLE 1. Summary of the systematic errors 

I Source I RUN-I(ppm) I RUN-II(ppm) 



TDC non linearity 


< ± 15 


< ± 15 


Pile up 


< + 10 


<+ 10 


Pickoff estimation 






(l)3y subtraction 


< ± 89 


<±91 


(2) Ge efficiency 


< ± 33 


<±28 


(3) YAP efficiency 


<±64 


< ± 19 


Physics source 






( 1 ) Zeeman effect 


< - 5 


< - 5 


(2) Stark effect 


< + 3 


<+4 


(3) 3y annihilation 


< - 91 


< - 33 



Total -147 and +115 -104 and +98 



physics source(Stark effect and 3y annihilation in the pick-off collision). 

The predominant contribution to total systematic error is produced by uncertain nor- 
malization of the 3y spectrum. That is, the number of pick-off events are determined 
by subtracting the normalized 3-y spectrum of Monte Carlo simulation from the o-Ps 
spectrum, where changing the normalization factor affects the ^ P i c k{t)/My values and 
eventually propagates to the final result. Since the sharp fall-off of the 3y-spectrum at 
511 keV is solely produced by the good Ge energy resolution of o = 0.5 keV, this 
subtraction only affects the lower side of the pick-off spectrum such that improper sub- 
traction results in asymmetry of the pick-off spectrum shape. Comparison of the asym- 
metries of the pick-off peak shape and the prompt peak annihilation spectrum is a good 
parameter for estimating this systematic error. The errors are about 90 ppm for both 
measurements. 



The Stark shift stretches the lifetime of Ps atoms, i.e., a probative calculation shows 
that the shift is proportional to a square of the effective electric field E such that 
AX-iy/X-iy = 248 • {E/E ) 2 , where E = m 2 e e 5 /h 4 « 5.14 x 10 9 V/cm. Silanol functional 
groups on the surface of the grain behave as an electrical dipole moment creating 
an effective field around the grains. Average densities of Silanol are measured to be 
OAAnmT 2 and the effective field can be analytically calculated such that the contribution 
to the o-Ps decay rate is determined to be —5ppm. These estimations were confirmed 
using results from precise hyper-fine-structure (HFS) interval measurements of ground 
state Ps in silica powder lllOll . where the interval is proportional to the size of Stark effect. 
Considering the difference in powder densities used, the HFS results are consistent with 
our estimation. 

Other sources of systematic errors: Error contribution due to the Zeeman effect is es- 
timated using the measured absolute magnetic field around the positronium assembly 
(—5 ppm). Since the 3-y pick-off process can only occur at a certain ratio, the calculated 
relative frequency Osyjoj.y ~ 1 /378 is consistent with previous measurements flUl . be- 
ing —91 ppm for RUN-I and —33 ppm for RUN-II. 

The above discussed systematic errors are regarded as independent contributions such 
that the total systematic error can be calculated as their quadratic sum, resulting in 
— 147 ppm, +115 ppm for RUN-I and — \04ppm, +9Sppm for RUN-II. The combined 
result with systematic error is X-iy = 7.0401 ± 0. 0006 (st at .^qq^ 1 (sys .) /is -1 , and total 
error is 150 ppm. 



Results and discussion 

Figure 5 shows the summary of the measured decay rate after 1995, three results are 
obtained with our method and one result is obtained using thin polymer, with which the 
produced positronium is very slows and the the effect of the thermalized positronium is 
suppressed. These four results are consistent with each other, and the combined value is 
X?,y = 7.0401 ±0.0007 /is -1 and shown in the red arrow in the figure. The correlations of 
the systematic errors are carefully taken into account. The combined result is consistent 
with the all four results and accuracy is 100 ppm. This result is consistent with the O(oc 2 ) 
correction[? ], and differs from the only up to O(a) by 2.6a. 

HYPER FINE SPLITTING 
Current status 

The hyper fine splitting (HFS) is the difference of the energy level between the singlet 
(l 1 ^) state, parapositronium(p-Ps) and the o-Ps. The HFS of the positronium is large 
because of the following two reasons: (1) The magnetic moment is proportional to the 
inverse of the mass, and large spin-spin interaction is expected for o-Ps. (2) o-Ps has the 
same quantum number as photon, and o-Ps has quantum oscillation into virtual photon, 
o-Ps — » 7* — > o-Ps. This frequency of 87GHz contributes only to the o-Ps, and makes 
HFS larger. 
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FIGURE 5. Decay rates measured in the last four experiments are listed chronically. The last arrow 
shows the combined result of the last four measurements. The red and blue lines shows the QED prediction 
calculated up to 0(a 2 ) and O(a), respectively 



The precise measurement of the HFS gives the direct information about the higher 
order calculation of the QED, especially of the bound state QED. If an unknown light 
particle (like as axion or millichaged particle) exits, it contributes to the energy level, 
and makes discrepancy from the QED prediction. Since the quantum oscillation has 
good sensitivity to the hypothesis particle, whose coupling is super-weak, the precise 
measurement of HFS is good tool to search for the new physics beyond the Standard 
Model indirectly. 

The precise measurements have been performed in 70's and 80's. These results are 
consistent with each other and the final precision is 3.6 ppm tl2D. No higher order 
correction(> 0(a 2 )) has been performed at that time and these results "were" consistent 
with the 1st order calculation of the QED. New method to calculate the higher order 
correction on the bound state is established in 2000, and the 2nd and 3rd corrections 
have been performed[13]. The QED prediction is 203.3917(6) GHz and it differs from 
the measured value of 203.38910(74) GHz. The discrepancy of 3.5rj is observed. It is 
still marginal for statistical fluctuation, but there would be possibility of the new physics 
or the common systematic errors in the previous measurements. (There is also possibility 
that NRQED approximation is wrong.) 



Old experiments and the systematic errors 



In the previous all measurements, the HFS transition was not directly measured, since 
203GHz was too high frequency to be handled. The static magnetic field makes Zeeman 
mixing between m z =0 state of o-Ps and p-Ps, the resultant energy level of m z =0 state 
becomes higher than m z = ±1. This energy shift is proportional to the HFS energy level. 
Figure 6 shows the energy levels as function of the static magnetic field. The m z = ±1 
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FIGURE 6. The energy level of Ps as function of the static magnetic field. Av represents the HFS, and 
o-Ps state ( 3 Si) is higher than p-Ps (^o). Arrow shows the Zeeman shift between M- = and M z = ±1 at 
B=8kGauss. 



state of o-Ps do not couple to the static magnetic field, and do not change the energy 
level. 

Positronium was produced with the /3 + source and gas(/V2, Argon etc) in the RF 
cavity, in which the high power 2.3GHz microwave was stored. Static magnetic field 
(about 8kGauss shown with arrow in the figure) was also applied and scanned near 
the resonance of the Zeeman shift. When the Zeeman energy shift is just on the RF 
frequency, the transition from m z = ±1 to m z =0 increases, and the m z =0 state decays 
into 2y immediately through p-Ps state. The 2y decay, which were tagged with back- 
to-back topology, increases on the resonance. HFS can be determined as center value 
of the resonance peak by scanning the magnetic field. Since Ps is produced in the gas, 
the Ps collides with the gas molecule and the electric field of the gas molecule makes 
the shift of the energy state, which is called as the Stark effect, and about lOppm shift 
was observed for 1 atm gas [|12|l . The HFS in gas were measured by changing the gas 
pressures. The measured values were extrapolate to zero pressure, and the HFS in the 
vacuum was obtained. This extrapolation method was the exactly same as in the decay 
rate measurements |Q1 0, [H] . 

There are two possibilities of systematic errors in this method. (1) First is the effect 
of the non-thermalized o-Ps, which is the same as in the decay rate measurement. 
Extrapolation procedure is assumed that Ps is well thermalized and the mean velocity of 
the Ps is the same for the various pressured gases. We have already shown in the decay 
rate measurements that this well-thermalized assumption is not satisfied, and that this 
method makes big systematic errors. The non-thermal Ps would affect also on the HFS 
measurement. (2) The Ps is widely spread in the RF cavity (size 17cm in diameter) and 
non- uniformity of the static magnetic field was the main source of the uncertainties of 
the results. Size of the used magnet were limited and the uniformity of the filed was 
about lOppm level, which has been corrected. 



FIGURE 7. Superconducting magnet: Bore size is 80 cm and the uniform magnetic field can be applied 
up to IT. 



Our new methods 

We propose the following two new methods in order to solve these systematic uncer- 
tainties in the previous measurements. 



Conventional Zeeman method using the magnet field 

22 Na positron source(lMBq) is installed in the thin plastic scintiHator(100/im), and 
the timing information of the positron emission is tagged (positron emits at t=0). There 
are two benefits to have time information: (1) Many part of positron is just annihilated 
into 2y at t=0, and this 2 7 annihilation events make S/N seriously worse. The 2y annihila- 
tion background events can be removed dramatically with the requirement of t>10nsec. 
(2) Thermalization process can be determined with the energy spectrum measured by 
the germanium semiconducting detectors, as the same as in the measurement of the de- 
cay rate. Stark effect (material effect) can be measured changing gas pressure and the 
non-thermalization effect are measured directly. 

Magnetic filed is provided with the superconducting magnet developed for the medi- 
cal NMR. Figure 7 shows the photograph of the magnet, which has the large bore size 
(80cm) and strong (about IT) magnetic field. The magnetic filed is uniform well, we 
can reduce the uncertainties due to non-uniformity of the magnetic filed. 0.9T of the 
magnetic filed will be applied with 2.8GHz RF (power 500W) microwave. Cavity with 
TM1 10 mode resonance is used to store the RF power. 

Figure 8(a) shows the gamma ray detectors and the RF cavity. Six LaBr3 scintillators 
and two germanium detectors are used to detect gamma ray from the Ps decay. LaBr3 
scintillator has the fast signal shape as shown in Fig. 8(b). Fast raising time, faster than 
10 nsec, provides the good time resolution of about 300psec (FWHM at 51 IKeV), and 




FIGURE 8. (a) The experimental setup of the gamma ray detectors: 6 LaBr3 scintillators and 2 Ge 
detectors are located around 2.8GHz RF cavity. Since these are located in the magnetic filed, the fine- 
mesh PMTs are used and the angel between the magnetic field and the PMT axis is set at 30°, at which 
the magnetic field effect on the PMT gain is about 10. (b) The observed signal shape of the LaBr3 for 
667KeV gamma ray emitted from 137 Cs. 



TABLE 2. The expected numbers of events for 2y 
and 3 7 decay with the both tagging method. Numbers 
will be detected on the Zeeman resonance for one day 
run 





geometrical tagging 


Energy tagging 


2y decay 


4.3 xlO 4 


3.2xl0 6 


37 decay 


7xl0 3 


2xl0 6 



the fast decay time (26nsec) grantees the pileup effect can be reduced. Energy can be 
measured with the high resolution (FWHM=3% at 511keV) as shown in Fig. 8(b). The 
energy resolution is better than Nal by factor 2, then 2y processes can be tagged with 
the single scintillator using the energy information (Energy tagging). Since the tagging 
efficiency of the energy tagging is significantly higher than the 'geometrical tagging', 
in which 2y decay is tagged with two scintillators located back-to-back. The expected 
numbers are summarized in table 2. First measurement with this new method will be 
performed in this Autumn, and a few ppm accuracy will be obtained for a half year. 



Direct transition using sub THz RF 

(Sub) THz RF is the unexplored filed, since (sub) THz RF has the characters between 
the optics and the radio wave. Developing the high quality/high power source for the 
(sub) THz RF is useful and interesting for both science and technology. We can observe 
the direct transition of the HFS without magnet field, if we have good source of 203 GHz 



FIGURE 9. (a) Photograph of Gyrotron FU CW-4: High power (500W) and high frequency (400GHz) 
sub THz wave can be produced, (b) The expected resonance curve with Q=6000 cavity and 100W RF 
source. 



RF. It is free from the systematic error due to magnetic filed, and it open the new 
paradigm of the atomic/particle physics using (sub) THz. 

Figure9(a) shows photograph of the gyrotron RF generator[14], in which high 
power(>100W) sub THz RF can be produced continuously. The accelerated electron 
emits RF in the strong magnetic filed, and the one/specific frequency is enhanced with 
the resonance cavity. Arbitrary RF frequency can be picked up with the tuning reso- 
nance cavity. Tunable frequency is also big challenge of THz RF source. Improvement 
stability of the power is in progress and the final design of stability is about 10 ppm. 
This is big challenge on the THz technology and we can expect various spin off for not 
only basic science but also technology. 

203 ± 1GHz (tunable frequency) CW RF source will be ready in this year. The gamma 
ray detectors mentioned above are used, and the positron timing is also measured in 
order to reduce annihilation 2y and measure the thermalization process. The Fabry- 
Perot cavity is used for RF storage, since 203GHz has optical characters. Figure 9(b) 
shows the expected resonance curve with 100W RF power and Q=6000 cavity. The 
width(FWHM 0.63GHz) is determined with the decay rate of p-Ps, and the center value 
can be determined with accuracy of about 100 ppm. 

CONCLUSION 

The decay rate of o-Ps was measured using a direct 2y correction method in which 
the thermalization effect of o-Ps is accounted for and integrated into the time spectrum 
fitting procedure. The new result is A 3y = 7.0401 ±0.0006(stat.)±Q$j^(sys.) ns~ l , and 
total error is 150 ppm[9]. It is most accurate result and agrees well with the last three 
results ASS]. The combined result of these four measurements is 7. 0401 ±0.0007 /is , 
which is consistent well with th NRQED prediction corrected up to 0(a 2 ) term|4|], and 



differ from the result up to 0(a) by 2.6a. It is the first result to validate the 0(a 2 ) 
correction. 

The Hyper Fine Splitting between o-Ps and p-Ps is sensitive to the higher order 
corrections of the QED calculation, and the discrepancy of 3.5c is found between the 
measured values and the 3rd order QED prediction. We point out that there are two 
possibilities of the common systematic errors in the previous all measurements. (1) Non- 
thermalized Ps makes error in the extrapolation method as the same as in the decay 
rate measurements. This was serious source of the orthopositronium lifetime puzzle. (2) 
Non- uniformity of the magnetic field is also serious source of the systematic error. 

We propose new methods to measure HFS precisely without the these uncertainties. 
Thermalization processes are measured with the 2y/3y ratio, and also 2 7 annihilation 
background events can be reduced by requirement of t > 10 nsec. LaBr3 scintillators are 
used for tagging 2y decay, since its energy resolution is better than the conventional scin- 
tillator i.e. Nal and Csl, we can tag 2y decay only with the energy information, not with 
back-to-back topology. It gains by factor 70 for the tagging efficiency. Large bore su- 
perconducting magnet is prepared for new measurement, it provides the uniform/strong 
magnetic filed. 

203GHz high power and stable RF source are developed. It makes possible to transit 
p-Ps into o-Ps directly. We can measure HFS directly instead of Zeeman shift. It will be 
the first application of sub THz RF source to the basic science. 
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